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^vq ■ effective quark-meson model. We investigate the transition between infinite-volume 



We study the scaling behavior of the two-flavor chiral phase transition using an 



and finite- volume scaling behavior when the system is placed in a finite box. We can 
estimate effects that the finite volume and the explicit symmetry breaking by the 



Q^' current quark masses have on the scaling behavior which is observed in full QCD 

cd : 

lattice simulations. The model allows us to explore large quark masses as well as the 
chiral limit in a wide range of volumes, and extract information about the scaling 

> 

I/-) ' regimes. In particular, we find large scaling deviations for physical pion masses 



and significant finite volume effects for pion masses that are used in current lattice 
simulations. 

PACS numbers: 12.38.Gc, 64.60.ae, 64.60.an 



^ ■ I. INTRODUCTION 



The order and the exact nature of the chiral phase transitions in QCD is important for 
our understanding of the dynamics at the early stages of the universe and of heavy-ion 



collisions, see e. g. 



2] . However, the order and nature of this phase transition has proven 



to be notoriously difficult to pin down |3jjl4(. It has long been surmised, from symmetry 



15], that the phase transition in QCD at 



arguments and a renormalization group analysis 
physical quark masses is a crossover and in the domain of a critical fixed point of a theory 
with 0(4) symmetry in d — 3 dimensions. In the limit of Nf = 2 massless quark flavors, a 
second-order phase transition governed by this fixed point is expected, whereas the phase 
transition becomes first order for Nf = 3 massless flavors. 
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Recently significant progress has been made towards establishing the scaling behavior for 



;he two-flavor transition 



8, llj and the 2 + 1-flavor transition with one heavy quark flavor 
121 ] . The transitions for both theories are expected to fall into the 0(4) universality class. 

Since the chiral flavor symmetry is broken explicitly by the quark masses, it is difficult 
to directly observe the critical behavior in lattice simulations. The transition becomes a 
crossover and the behavior is only pseudo-critical. The actual critical behavior can therefore 
only be ascertained by means of a scaling analysis. 

Because of the similar numerical values of the critical exponents for 0(N) models, see e. g. 



16 



26], it is also difficult to distinguish the expected 0(4) scaling behavior from 0(2) scaling 



behavior, which might apply due to a residual symmetry in the staggered implementation of 
fermions in lattice QCD. In this context it is useful to also compare results from simulations 
to scaling functions from the appropriate models 



Scaling in the 0(4)-model 27H32] and in quark-meson models with 0(4) symmetry 



33 



361 ] in infinite volume have been investigated previously using various methods. These studies 
provide results for the scaling functions in infinite volume, which can be used for comparisons. 
Most of these investigations have focussed on the scaling regime close to the critical point, 
with extremely small amounts of explicit symmetry breaking, and have not investigated how 
this scaling behavior is connected to the behavior at the physical pion mass. In order to 
assess the feasibility of a scaling analysis on the lattice at physical values of the pion mass, 
it appears useful to perform such an investigation in the context of a model for dynamical 
chiral symmetry breaking, e. g. NJL-type models. In such models physical pion masses and 
the chiral limit are both readily accessible. 

Universality arguments for the scaling functions do not, in general, allow us to make 
observations away from the critical region where the universal behavior obtains. They are 
only applicable where the physics of the systems are dominated by the critical long-range 
correlations. Probing shorter distance scales away from the critical point, aspects of the 
underlying short-range physics once again emerge, and different systems once again behave 
differently, even when they are in the same universality class concerning the critical behavior. 
However, the size of the 'critical regime' depends on the additional scales which come into 
play in different theories. For example, in NJL-type models as well as in QCD, there is 
a gap between the pions as the light Goldstone-modes, and the next heavier excitation. 
This makes e.g. chiral perturbation theory viable as a low-energy effective theory of QCD. 
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For this reason, we expect that for values of the pion mass below the physical values and 
temperatures much smaller than the critical temperature, differences between the behavior of 
the models and QCD ought to remain small, so that the model results still have relevance for 
our expectations about full QCD. In the vicinity of the critical temperature this assumption 
may no longer hold since at the transition lots of higher resonances come into play j^ . 
In fact, the dependence of the critical temperature on the pion mass is much stronger in 
NJL-type models than that which has been found in lattice simulations. We shall address 
this issue below. 

In addition to the issue of finite pion masses, lattice simulations are also always performed 
in finite simulation volumes. Thus one has to take care to exclude finite-volume effects in 
order to observe the scaling behavior expected for infinite-volume systems. Such finite- 



volume effects can also affect the phase transition at finite density 



38|-|40|. Alternatively, 



one can turn the appearance of finite-size effects into an advantage and perform a finite- 



ly y, H : 



size scaling analysis 4lL |42| to obtain additional information about the scaling behavior 



43 



445| . Finite-Size scaling functions relevant for the analysis of the chiral order 



parameter have been obtained e.g. in 44, 



45]. 



We can distinguish between two regimes, one where the scaling behavior conforms to that 
for infinite volume, and one where one observes clear finite-size scaling effects. In between 
these two regimes, it is not obvious how the finite volume will affect the scaling behavior, and 
an investigation of the transition between these regimes in a model system is worthwhile. 

A crucial part of the underlying argument for our expectations of the scaling behavior at 
the chiral phase transition in QCD relies on the idea of dimensional reduction: If we consider 
a system in infinite volume, but at finite temperature, as a system in Euclidean space, close 
to a critical point, then the wavelength of the critical fluctuations eventually become larger 
than the extent of the system in the Euclidean time direction. This leads to an effectively 
three-dimensional system, and subsequently to our expectation of three-dimensional scaling 
behavior in the 0(4) universality class for two-flavor QCD. 

In order to observe infinite-volume scaling behavior in the expected universality class in 
finite- volume lattice simulations, on the one hand the extent in Euclidean time direction has 
to be small enough to lead to dimensional reduction. On the other hand the extent in the 
spatial directions has to be large enough to minimize finite-volume effects. The scale set by 
the temperature must then be compared to that set by the long-range correlations, which is 
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bounded by the pion mass. This is a potential issue for lattice simulations, where the overall 
scale for the lattice spacing is set by the value of the gauge coupling, but the aspect ratio 
of the lattice in spatial and Euclidean time direction remains fixed. It appears therefore 
worthwhile to investigate both the infinite-volume scaling behavior and the deviations from 
this behavior for fixed aspect ratio in a finite volume, both for very small pion masses close 
to the chiral limit, and for physical values. 

Considering that the physical pion mass is around 140 MeV, and the corresponding chiral 
phase transition temperature is approximately 150 — 180 MeV @, [hJ Q, Q > the proposition 
that the system should experience dimensional reduction for a physical choice of parameters 
does not appear a priori obvious. 

The scaling analysis presented in this work may help to shed light on the scaling analysis 
in actual QCD studies since the physics of the chiral phase transition is mainly determined by 
the long-range effective degrees of freedom at low momentum scales, namely the pions as the 
Goldstone modes of chiral symmetry breaking. As far as these degrees of freedom determine 
the behavior of the transition, according to universality arguments, the much simpler model 
system should exhibit the same critical behavior as QCD. Of course, our model approach 
cannot answer questions outside the applicability of the model. For example, the order of 
the phase transition is in our case already fixed by the 0(4)-symmetry of the model, while 
the order of the transition in QCD has not yet been unambiguously determined 
Moreover, we must also limit our investigation to the chiral phase transition. On the other 
hand, in calculations based on renormalization group methods the model has been used 
to investigate chiral properties of QCD beyond the mean-field approximation, such as the 



critical behavior and the quark mass dependence of the chiral transition 33|, |34|, |46( as well 



as the critical behavior at finite density 



471 ] . In the past few years it has also been combined 



with Polyakov loop results from lattice QCD simulations to improve the description of 
thermodynamical observables, see e. g. Refs. 0B|- 

Corrections beyond the mean-field 



approximation in the (P)NJL/(P)QM model have been considered in 551458]. 

The paper is organized as follows: In Sect. [Til we briefly discuss the model and our RG 
approach to a scaling analyis in infinite and finite volumes. In Sect. II III we present the 
results from our scaling analysis in infinite volume while we discuss scaling in finite volumes 
in Sect. IIVI Our conclusions are summarized in Sect. [V] 



5 



II. CHIRAL MODEL AND NON-PERTURBATIVE RG APPROACH 

In this section we discuss our RG approach to a finite- volume scaling study of the chiral 
phase transition in QCD. We briefly introduce the (linear) quark-meson model and discuss 
the derivation of the flow equations for finite and infinite volume studies. A detailed discus- 
sion of the derivation and the approximations involved can be found in Refs. 3l|, |46|, |59|, |60 |. 

To study the chiral phase transition and its scaling behavior for infinite and finite volume, 
we employ the chiral quark-meson model. This model is an 0(4)-invariant linear u-model 
with Nf = 2 quark flavors with N c = 3 colors and Nj = 4 mesonic degrees of freedom. The 
mesons are coupled to the (constituent) quarks in an SU(2)l x SU(2)r invariant way. We 
stress that it is an effective low-energy model for dynamical spontaneous chiral symmetry 
breaking at intermediate scales. However, it does not contain gluonic degrees of freedom 
and is not confining. 

At the scale A, the quark-meson model is defined by the bare effective action 

r A [q,q,4>] = j d^^ + g^a + if-n^^ + ^d^f + U^-Ha] (1) 

with <f) T = (a, 7?) and ^ denote the fermion spinors associated with the quark fields. 
We choose the first component of the vector <p to be the radial mode associated with the a 
meson. Note that due to the explicit symmetry breaking the ground state of the theory is 
only symmetric under 0(3) transformations. The mesonic potential at the ultraviolet (UV) 
scale is characterized by two couplings, m\ and Aa, 

U^)= l -ml^ + \\ K ^)\ (2) 

The quarks and the mesons are coupled via chirally symmetric Yukawa term with g being 
the coupling. The linear term in a results from a bosonization of the current quark mass 
term ~ fm c \l/ and the symmetry breaking parameter H is therefore related to the current 
quark mass and m(A) = m\: H = m c m\/ g. We study the quark-meson model in the so- 
called local potential approximation (LPA), where we neglect a possible space dependence 
of the expectation value ((f)) and take the wave-function renormalizations and to be 
constant, Z^ = 1 and Z$ — 1. This approximation should not be confused with a mean- field 
approximation. In fact, our approximation includes already beyond mean- field effects. A 
detailed discussion of the relation of the present approximation (LPA) to the mean-field 
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approximation can be found in Refs. 6l|, |62|. At finite temperature we moreover neglect 



a possible difference of the wave-function renormalization parallel and perpendicular to the 
heat-bath. It has been found that the latter approximation does not strongly affect the 



dynamics near the phase transition 62[|. Since the anoma 



and are small compared to one, see e. g. Refs. [33 



ous dimensions associated with 



62 



63], our approximation, in 



which the running of the wave-function renormalization is neglected, is justified for a study 
of finite-size scaling. In fact, corrections beyond the local potential approximation changes 



the resulting critical exponents only at the one-percent level [26j, [33|, |63j . For products of 
critical exponents that enter our scaling analysis the changes are slightly smaller. 

For our derivation of the RG flow equation for the effective action we employ the Wetterich 
equation [64j |: 



d t T k 



I 



STr {d t R k ) 



,(2) 



(3) 



where the dimensionless flow variable t is given by t = lnf fc/A). R eviews of and introductions 



to functional RG approaches can be found, e. g., in Refs. 



65 



- l74j . The regulator function 



specifies the details of the Wilsonian momentum-shell integrations and has to satisfy certain 



constraints 



64j. Since the choice of the regulator function is at our disposal, we can use it 



70 



75 



to optimize the RG flow 

Rk(p 2 ) = p 2 r(p 2 /k 2 



77l | . In the following, we employ 77 1 



with 



r(x) 



1 6(1 



x) 



(4) 



In order to derive the RG flow equations for a system in a finite four-dimensional Euclidean 
volume L 3 x 1/T with temperature T, we replace the continuous momenta by discrete 
momenta and correspondingly each momentum integral in the evaluation of the trace in 
Eq. ([3]) by a sum: 

3 oo oo oo 



p -> 47T n = Air (nf + n 2 + n 



and 



d 3 p 



t) e E E • 

n±=— oo ri2=— oo n:j=— oo 



Since we are ultimately interested in a comparison to scaling behavior of the chiral order 
parameter as measured in lattice QCD simulations, see e. g. Ref. 12[, we choose periodic 
boundary conditions for the bosons and fermions in the spatial directions. 

For studying scaling behavior it is convenient to deal with dimensionless quantities rather 
than dimensionful quantities. Therefore we introduce the dimensionless potential u, the 
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dimensionless fields (p and the dimensionless symmetry breaking parameter c by 

Uk = k~ 4 Uk , <fi = , and c = k~ 3 H. (6) 

Note that the Yukawa coupling in d = 4 is already a dimensionless quantity. With these 
definitions the flow equation for the effective potential of the quark-meson model in a finite 
box with length L at finite temperature T is then given by 



dtu k {^ 2 ) = -Au k + B{kL) 



2N c N f 



(l - n F (e q , fj) - n F (e q , -//) j 



(7) 



where the first two terms correspond to contributions of the mesonic modes, and the last 
term with opposite overall sign corresponds to the quark contributions. The effective energies 
are given by 



6i = y 1 + m? , ie {n, a, q} , (8) 

with 

2 # u fc 2 o^Wfc 2 <9 2 U fc 2 2 2 /„n 

^=V' mCT = V +4 ^w' m «=^- (9) 

The temperature dependence of the RG flow of the effective potential is governed by the 
bosonic and fermionic distribution functions ub and np, respectively: 

nB ® = -^r—[> Wj(e '^ %fc-/o/f + i ' (10) 

where t = T/k denotes the dimensionless temperature. The dependence on the finite spatial 
volumes is encoded in the function B: 

1 oo oo oo 

W = 4 E E E e((*L) 2 -47r 2 n 2 ). (11) 

^ ' ni=— oo ri2=— oo n3=— oo 

The asymptotic behavior of this function for small and large dimensionless box sizes kL is 
given by 

lim B(kL) ~ — — |— r and lim B(kL) = At . (12) 
kL^o v ' (kL) 3 kL^oo y ' 6tt 2 v 7 



s 



The behavior for small kL reflects the fact that the dynamics of the system is mainly 
governed by the spatial zero modes in thi s limit. On the other hand we recover the same 
flow equation as found in Refs. 



47 



78|] for kL ->> oo. 

In order to solve the RG flow for the scale-dependent effective mesonic potential Uk, we 
expand the potential in a Taylor series in scale dependent local n-point couplings a n ^ around 
its scale dependent minimum (</? ) 



u k (ip 2 



n=0 



2\n 



(13) 



The presence of the symmetry-breaking term —Ha in our ansatz ([I]) induces a shift of the 



minimum from its value in the chiral limit. Following Ref. 





45|, the condition 



dif { 



M<P 2 ) = c (14) 

<PO=( t Po), t Pi=° 

ensures that the potential is always expanded around the actual physical minimum. From 
Eq. (I14p it follows that the RG flow of the coupling a^k and the minimum (tpo) are related 
by the simple condition 

Oi,fc(^o) = c. (15) 

This condition keeps the potential minimum at tp = ((</?o)>0). 

The RG flow equations for the couplings a n ^ and (tpo) can be obtained by expanding 
the equation for the effective potential, Eq. ([?]), around the scale-dependent (ip ) and then 
projecting it onto the derivative with respect to k of the ansatz, Eq. ([TBI . This procedure 
results in an infinite set of flow equations for (tpo) and the couplings a n (k). In order to 
obtain a finite set of flow equations, we truncate the Taylor series, Eq. ([TBI , at a fixed order 
-^Vmax = 4 and include thus fluctuations around the physical ground-state configuration up to 
order 2iV max in the fields. Note that such an expansion represents a systematic expansion in 
m-point functions r^ m ^ where m = 2n determines the number of external legs. The quality 
of such an expansion of the order-parameter potential in powers of ip 2 has been studied 



quantitatively in Ref. 



63] at vanishing temperature and for the proper-time RG in LPA at 



finite temperature in Ref. 79|. Moreover the order-parameter potential has been computed 
in LPA without making use of a Taylor expansion in ip in Ref. 
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III. SCALING ANALYSIS IN THE INFINITE VOLUME LIMIT 

In the vicinity of a critical point, where the dynamics of the system are dominated by 
critical long-range fluctuations, the singular part of the free energy density of the system 
satisfies to leading order the scaling relation 

f s (t,h)=£- d f s (tiy\h£^), (16) 

where £ is a dimensionless rescaling factor which can be chosen arbitrarily, and t = (T — 
T c )/T and h = H/H are the reduced temperature, measured from its critical value, and 
the external symmetry-breaking field, normalized in a suitable way. 

The exponents y t and y t specify all critical exponents for the scaling behavior, 

1 P S 

yt = -, Vh = — , (17) 

V V 

when taken in combination with the additional scaling relations 7 = (3(5 — 1), 7 = {2 — rf)v. 

As a consequence of the scaling relation, observables such as the order parameter, iden- 
tified with the pion decay constant M = f n in the model, and the susceptibilities Xn f° r 
transverse Goldstone modes and \ G for longitudinal modes can be expressed in terms of 
universal scaling functions. By choosing the scaling factor £ such that either t£ yt = 1 or 
hl Vh = 1, the free energy density becomes a function of only a single scaling variable, with 
an explicit dependence on either t or h. Thermodynamic observables which can be expressed 
in terms of derivatives of f(t, h) with respect to its arguments can then also be expressed in 



terms of such scaling functions 80]. 



For the order parameter, one finds the scaling relation 

M = h l ' s f M {z) 1 z = t/h 1 ^ (18) 

where z is the scaling variable, and /m(^) is the scaling function normalized to /m(0) = 1. 
Asymptotically for small values of h and t < 0, i.e. for large values for —z, we have 
/^(z) ~ (— z) 13 . These two conditions determine the normalization constants To and Ho, 
such that M = h l / s for t = and M = (-tf for h = and t < 0. 

In the LPA used in this paper, the static susceptibilities are related to the masses of the 
mesonic modes according to 
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for the transverse and the radial modes, respectively. Since the susceptibility for the trans 
verse mode is related to the order parameter M = (<p ) = f n according to Eq. ( 1T5|) . i.e. 

M (<p ) 



(20) 



H H ' 

the transverse susceptibility does not contain any additional information beyond that con- 
tained in the behavior of the order parameter, and we will therefore not consider it separately. 

The longitudinal susceptibility Xa — fff can be expressed in terms of the scaling function 
f x (z), which is related to the scaling function Jm{ z ) an d its derivative according to 



Xa 



■fxW 



Mz) - jf M {z) 



(21) 



H Jxy ' H 5 

It corresponds to the chiral susceptibility, i.e. to the susceptibility of the chiral condensate 
(i/jiI') with respect to a change in the current quark mass m c . 

We wish to stress that the scaling analysis for a theory with d = 3 is here performed 
for a theory in d = 4 Euclidean dimensions, where the temperature (in a field-theoretical 
sense) is given by the Euclidean time extent of the volume. This means that we will in 
fact only observe scaling in a three-dimensional universality class when the conditions for 
a dimensional reduction are met. In contrast, many earlier determinations used a three- 
dimensional theory, where the transition was determined bv the critical value of one of the 



couplings of the theory see e.g. [3l|, |35| for RG and [28 



44 



8 If for spin model results. A 



scaling analysis in infinite volume based on functional RG approaches, in which temperature 



33 



35 



36]. 



has been introduced in a field-theoretical sense, has been performed in Refs. 
In this work, we study scaling for small and large pion masses in infinite but also in finite 
volumes. 

In the LPA, the Yukawa coupling g, the symmetry breaking field H as well as the expan- 
sion coefficients a n & of the order parameter potential Uk are parameters at the UV cutoff 



scale which are at our disposal. In principle it is possible to fix t 
principles by employing an RG group approach to full QCD [61 



leseparameters from first 



82 



86j. However, we shall 



not follow this strategy in the present paper. Since we are interested in a study of the chiral 
phase transition in QCD, we fix the parameters such that we reproduce the physical values 
of the pion decay constant, the pion mass and the constituent quark mass in the infrared 
(IR) limit, i. e. for k ->• 0: 



M w w 138 MeV, / w w 93 MeV, 



M q w 



298 MeV. 



(22) 
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Since we use only three IR observables to fix at least three UV parameters, there is some 
ambiguity in the parameter-fixing procedure. However, we have checked that our results 
depend only weakly on the actual values of the UV parameters of our model, provided a 
given parameter set yields the same IR values for our physical observables. This observation 



is in accordance with earlier studies, see e. g. Refs. 31 



451 ] . To be specific, we choose 



2 

= 3.2, «i,a = ^ ^ 0.547, fl2,A = 2Aa ~ 67.2 , a n>A = for n > 3 (23) 

at the scale A = 3 GeV. We have also checked that our results for, e. g., the critical temper- 
ature are independent of A for A > 3 GeV. For A < 3 GeV we find that our results depend 
on A. Our choice for A is larger than typical values for the UV cutoff in (P)NJL/(P)QM- 
type model studies where A is identified with some hadronic scale below which a mesonic 
description of QCD may become applicable. A (weak) dependence of, e. g., the critical 
temperature on A does not usually play a role in such model studies. In the present study, 
however, we have to take care that our results for, e. g. the critical temperature, do not 
exhibit a dependence on the UV cutoff since it would spoil our scaling analysis. From a 
phenomenological point of view, choosing a large UV cutoff A for a low-energy model is 
not an issue, provided the UV parameters of the model have been chosen such that the IR 
physics remains unchanged. The predictions from the model for, e. g., the (chiral) phase 
boundary, are then not (or only slightly) affected by the choice for the cutoff A. 

From now on we leave our choice for the couplings at the initial UV scale unchanged for 
all temperature and volume sizes. We are then left with one parameter, namely the external 
symmetry breaking field H, which mimics the current quark mass. Different values of H 
translate directly into different values of the pion mass M n . For H — > (chiral limit) we 
find the following value for the critical temperature: 

T c « 144.949346731961 MeV . (24) 

Note that in our model study such a high accuracy for T c turns out to be necessary to resolve 
the scaling region in the infinite- volume limit, see below. This has also been observed in 



scaling studies of O(N) models in d = 3, see e. g. [3l|, |45|. Studying the model for various 
values of H in the vicinity of H = 0, we obtain from Eq. (jl~8]) the values for the normalization 
constants To and Hq: 

T » 23.862066412513776 GeV, H « 346.37722832663883 GeV 3 . (25) 
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We find that the critical exponents of our model are consistent with the well-known LPA 



values found in studies of O(N) models [20|, |3l|, |45|: 



/?w 0.4022, 5^5.0. (26) 

In the following we shall use these values for the critical exponents in our scaling analysis. 

Let us now briefly discuss our results in the infinite- volume limit. In Fig. [T] we show 
the chiral susceptibility x<* and the rescaled chiral susceptibility for various pion masses 
from M n = 0.2, . . . ,0.9MeV. From the results for the chiral susceptibility Xa, (left panel) 
we deduce that the scaling region in our model is indeed very small. This requires a high 
accuracy in the determination of T c , the critical temperature in the chiral limit. Since the 
susceptibility \ a is proportional to the squared correlation length, we define the peak of the 
susceptibility Xa to be the pseudo-critical temperature T P (H) = T p (M vr ) which is associated 
with long-range correlations. For high temperatures T ^> T c the system is outside of the 
scaling region and the results for the susceptibility Xa fall onto a single line, Xa ~ 1/T 2 . 

After rescaling, the curves for Xa for M n = 0.2, . . . , 0.9 MeV fall onto a single line and 
are almost indistinguishab 



at z p w 1.3155, see also 



31 



e at the scale of the plot. The maxima of the curves are located 
, 4|| . Thus scaling corrections are bound to be small in this pion 
mass range. However, corrections to scaling become soon apparent for larger values of the 
pion mass, as we shall see below. 

In Fig. [2] we present our results for the rescaled chiral susceptibility Xa an d the rescaled 
pion decay constant as a function of the scaling variable z for M n = 0.5 , 75 , 138 , 200 MeV. 
The pion masses used to obtain these results include the physical value as well as the 
currently smallest value used in lattice simulations 12|, |87J . We observe that the curves for 
the rescaled susceptibility and the rescaled order parameter do not fall onto a single line. 
This ought to be the case if corrections to scaling were small in this pion mass regime. On 
the contrary, the rescaled susceptibility and the rescaled order parameter differ significantly 
from the scaling function obtained for small pion masses. On the other hand the results 
for the rescaled susceptibility and the rescaled order parameter for M n = 75 , 138 , 200 MeV 
fall almost on one curve for z > — 1. This appears to signal the proper scaling behavior 
for these pion masses but only if we disregard the results for M n < 1 MeV. This might 
be helpful information for a scaling analysis in lattice QCD simulations. However, we have 
to keep in mind that our results have been obtained from a low-energy model in which 



13 




Figure 1: Left panel: Chiral susceptibility \a as a function of the reduced temperature t for various 
pion masses M n = 0.2, . . . , 0.9 MeV (from top to bottom). Right panel: Rescaled chiral suscep- 
tibility Xa as a function of the scaling variable z for various pion masses M n = 0.2, . . . ,0.9 MeV. 



0.4 I 1 1 1 1 1 1 1 1 1 2 



x 




-4 -3 -2 -1 1 2 3 4 5 -4 -3 -2 -1 1 2 3 4 5 

z = t/hVP 1 z = t/h 1 "* 



Figure 2: Left panel: Rescaled cr-suszeptibility as a function of z for various pion masses and 
L — > oo. Right panel: Rescaled pion decay constant f n (order parameter) as a function of z for 
various pion masses and L — > oo. 



the non-universal normalization constants are different from those found in QCD lattice 



simulations 



12 



87]. 



Finally we would like to comment on the pion mass dependence of the critical temperature 
in our model which is representative for NJL/QM models. From the definition of the scaling 
variable z and the universal peak position z p m 1.3155 of the rescaled susceptibility, we can 
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estimate the dependence of the pseudo phase-transition temperature T p on the pion mass 
for small values of H, i. e. for small current quark masses. From Eq. f fT8|) we deduce 



T P (H) = T p (M 7r ) = T c + T z p ( — 



= T C + c w Mf 



(27) 



where c v is a numerical constant which depends on the normalization constants and T P (H — > 
0) = T c . Since (35 ~ 2, we expect an almost linear dependence of T p on M-k-_ This is in 



accordance with the findings from lattice QCD simulations, see e. g. Refs. |88|, [89[. However, 
the values for the non-universal normalization constants To and Hq differ significantly from 
those found in lattice simulations. As a consequence, the slope of T p (M 7r ) is about one order 



of magnitude larger in chiral low-energy models than in lattice simulations, see e. g. [46 1 
for a more quantitative analysis of the slope. Thus, the universal properties of our model 
in the vicinity of the phase transition may agree with the findings from lattice simulations, 
provided two-flavor QCD falls into the 0(4) universal class. Non-universal aspects such as 
the pion mass dependence of T p seem to be incompatible with results from full QCD. The 
discrepancy in the slope c n can be traced back to the parameter-fixing procedure in our 
model approach: We have fixed the Yukawa coupling g, m 2 and A at the UV cutoff scale to 
reproduce the physical values for M n and M q for a given value of H. For our studies 
with various pion masses we have then only varied H but have left the UV values of g, m 2 
and A unchanged. However, these parameters have their own dependence on H (i. e. on 
the current quark masses) which is determined by quark-gluon dynamics at high momentum 
scales. We would like to add that the unknown dependence of the model parameters on the 
current quark masses may also affect the predictions of (P)NJL- and (P)QM-type mode 



for the phase boundary at finite temperature and quark chemical potential 50 



51 



53 



This dependence cannot be computec 



accessible in RG studies of QCD 
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with (P)NJL- and (P)QM-type models but is in fact 
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4-185] . A determination of the (current) quark mass 
dependence of these parameters from QCD RG flows is beyond the scope of this work and 
deferred to future work. 



IV. SCALING ANALYSIS IN FINITE VOLUME 

For systems in a finite volume, the critical behavior is modified because of the presence of 
the system boundaries. The linear extent of the volume L appears as an additional relevant 
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coupling: An actual critical point only exists in the limit 1/L — > 0. Phase transitions and 
the associated singularities in the free energy appear only in this limit. However, even away 
from the critical point, the behavior of the system will still be controlled by the critical fixed 
point. 

The critical behavior is affected as soon as the correlation length is of the order of the 
extent of the volume, and consequently the scaling regime is characterized by the ratio of 
correlation length and volume size. 

The singular part of the free energy density satisfies in leading order a scaling relation 

f s (t, h, L) = r d f s (t£ y \ M y \ LT\ (28) 

which now contains the volume extent L as an additional coupling. By choosing the rescaling 
factor t such that one argument is kept constant, we can use this relation to describe the 
system in terms of two variables only. In the limit of large volumes, the scaling behavior 
converges against the infinite- volume result. For this reason it is advantageous to choose 
the infinite- volume scaling variable z as one of these two variables. In leading order, we can 
then express the behavior of the order parameter M = f n as 

M{t,h,L) = L- p l v Q M {z,hL p& l v ), (29) 

i.e. as a function of the scaling variables z and hL 135 ^, where the volume dependence is 
now parameterized in the second variable. In order to reproduce the infinite- volume scaling 
relation in the limit L — > oo, the finite-size scaling function must satisfy the relation 

lim Q M (z,x) = x x/5 f M {z) with x = hL^ 5/u , (30) 

x— >oo 

where Jm{z) is the infinite- volume scaling function from Eq. (118)) . 

A similar relation can be derived for the longitudinal susceptibility, where in leading order 

X*(t,h,l) = L^Q x (z,hL^), (31) 

and Q x {z,hL l3S ' v ) is the finite -size scaling function for the susceptibility. Because of the 
scaling relations between the critical exponents, j/v = (2 — rj) is exactly 2 for our re- 



sults. For very small values of the symmetry-breaking parameter h, t 
Q x (z, hL@ s / u ) becomes essentially constant as a function of hL 135 ^ (see 
the susceptibility behaves as Xcr ~ L 2 for very small volumes. 



re scaling function 



45j). Consequently 
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Figure 3: Pion decay constant f w (order parameter) for z = (left panel) and z = z p (right 
panel) as a function of the normalized symmetry breaking parameter h for various values of the 
box length L. 



44 



45|. As 



The finite- volume scaling behavior of O(N) models has been considered in 
in the case of the infinite-volume scaling analysis, in these investigations relevant couplings 
which controlled the transition took the role of a temperature. In contrast, in the present 
investigation temperature is implemented as the finite extent of the Euclidean time axis in 
a four-dimensional Euclidean volume. 

In O(N) models as well as the chiral quark model considered here, at finite temperature 
the longitudinal correlation length is always bounded by the transverse correlation length. 
For this reason the scaling region can be characterized by the dimensionless product of box 
size and pion mass, M n L. 

In Fig. [3] we present our results for the order parameter f n in the quark-meson model 
in the finite-size scaling region. We have chosen to present results at fixed values of the 
scaling variable z. The choice z = corresponds to the critical temperature T = T c , and 
z = z p = 1.3155 corresponds to the position of the maximum (peak) in the longitudinal 
susceptibility Xa, which are shown in the left and right panel of Fig. [3j 

In agreement with our expectations Q, Q]' we find strong finite- volume effects for small 
values of the symmetry breaking, and a convergence to the infinite-volume behavior for 
strong symmetry breaking in the raw (unsealed) results shown in Fig. |3j We calculated the 
order parameter as a function of the symmetry-breaking parameter h for different values of 
the volume size, from L = 2 fm to L = 30 fm. For large values of the symmetry breaking 
the correlation length is small compared to the volume size, and the results converge to the 
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L [fm] 


2 


4 


6 


8 


10 


20 


30 


M w [MeV] 


308 


139 


85 


60 


45 


19 


11 


M W L 


3.12 


2.82 


2.59 


2.43 


2.30 


1.94 


1.75 



Table I: Pion mass M n (T — > 0, L — > oo) corresponding to the value of the scaling variable hL^ 5 / u 
at the bend point of the scaling curve in Fig. [H and corresponding values of the dimensionless 
product M^L. The values given are for the results with z = 0. 

infinite-volume behavior for large h. For decreasing h, the correlation length increases, and 
depending on the volume size L the deviation from the infinite- volume behavior sets in when 
the order of the correlation length approaches the volume size. For smaller volume size L, 
this happens for larger values of h. 

In Fig. H] the results for the finite-size scaled order parameter L^I V M = L^l v f n are shown 
as a function of the finite-size scaling variable hL 135 ^, for both z = and z = z p (left and 
right panel). The rescaled results fall onto a single scaling curve and thus show the expected 
finite-size scaling behavior. We can distinguish two different regimes in the rescaled results, 
where the rescaled results follow different power laws. 

For small volume size compared to the correlation length, the behavior is dominated by 
the effects of the finite volume. This part of the scaling function corresponds to those parts 
of the curves in Fig. [3] that deviate from the infinite- volume behavior. The "bend", where 
the slope in the double- logarithmic plot changes, characterizes the region in £/L where 
deviations from the infinite-volume behavior become the dominant effects. 

To give some guidance about the finite-size scaling region, in Tab. U we list the pion 
masses which correspond to the value of the symmetry-breaking parameter at the position 
where the slope changes. For a range of volume sizes, we give the value M 7r (T — > 0, L — > oo) 
of the pion mass at zero temperature and in infinite volume for the value of h at the bend 
point. We list in this example the values at the critical temperature, i.e. for the choice 
z = 0. 

For large values of the finite-size scaling variable hL^ s ^ u , i.e. for small £/L, the curve is 
characterized by the infinite-volume behavior. Due to the asymptotic large-volume behavior 
Eq. (l30i) . the finite-size scaling function for the order parameter behaves as x 1 ^ 6 for large 
values of the scaling variable x = hL 135 ^. In the double-logarithmic plot in Fig. HI the 
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Figure 4: Finite-size scaled order parameter f n L^^ u for z = (left panel) and z = z p (right panel) 
as a function of the scaling variable hL^^l v for various values of the box length L. 

exponent appears as the slope for large values of the scaling variable. We have checked 
explicitly that the inverse slope agrees with the value for the critical exponent 5 from our 
analysis, and find very good agreement for the results with z — 0. 

In the results from the 0(4) model presented in 45|, the asymptotic behavior for large 
volumes follows a power law with an exponent close to 1/5 for both z = and z — Zp as 
expected, with a slightly larger exponent for z = z p than for z = 0. In contrast, in the 
present case the exponent and thus the slope of the curve in the double-logarithmic plot for 
z = z p is somewhat smaller than for z = 0. The agreement between the finite-size scaled 
results for different values of L is better for z = that for z 



z p , which was also observed 



mg 



45] 



for the lattice spin model in [44|. The difference in the behavior of the asymptotic sea 
function between the quark- meson model in the present paper and the 0(4) model in 
is likely due to quark effects, as we will discuss below. 

As discussed in 45|], the presence of the additional coupling L requires the determination 
of an additional non-universal normalization constant L in order to determine a truly 
universal scaling function. A possible normalization condition is to require £ = L t~ u for 
h — 0. A direct comparison between the different scaling functions would be possible only 
after such an additional normalization. 

We will now turn from the dedicated finite-size scaling analysis to the finite- volume effects 
that appear if one performs a conventional infinite-volume scaling analysis in this volume 
and pion mass region. In Figs. [5] and [6] we show results for the susceptibility in a finite 
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Figure 5: Rescaled cr-susceptibility for M w = 48MeV (left panel) and M w = 75MeV (right panel) 
as a function of z for various box lengths L. 

volume for fixed values of the symmetry breaking parameter h, i.e. for fixed values of the 
pion mass M n . As examples, we have chosen the values M n = 48 MeV, M n = 75 MeV, 
M n = 138 MeV, and M n = 200 MeV. The parameters for the model are chosen in such a 
way that these values for the pion mass are obtained in the limit L — > oo and T — > 0. 

We present the results for the rescaled susceptibility XaH^h 1 " 1 ^ as a function of the 
infinite- volume scaling variable z = t/h 1 ^ 6 , as one would do in the absence of finite- volume 
effects. The purpose of this exercise is to investigate the deviations from the infinite-volume 
scaling behavior due to such effects and to estimate the volume size where they become 
relevant. This volume size is of course strongly dependent on the pion mass. 

As discussed above, for small volume size the susceptibility scales as \ a ~ L 2 since the 
correlation length is bounded in a finite volume, £ < L. Whether this behavior can be 
observed in the presence of explicit symmetry breaking depends once again on the dimen- 
sionless product M n L. For large symmetry breaking (M n L S> 1), this scaling behavior will 
not be evident. 

Starting with a comparatively small pion mass of M n = 48 MeV (left panel of Fig. |5J), 
we find that the expected decrease in the susceptibility with the volume can be observed 
only between the smallest volumes with L = 2 fm and L = 4 fm. This conforms to a rough 
estimate from the bound that the mass of the Goldstone mode places on the correlation 
length: a pion mass M n = 48 MeV corresponds to a length scale of approximately 4 fm. 
For larger volumes, the susceptibility is still bounded by the pion mass, i.e. the amount of 
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explicit symmetry breaking. The volume size becomes the limit only below this size. 

Nevertheless, we observe significant finite-volume effects in the susceptibility for larger 
volumes and pion masses. From Figs. [5] and El we see that for a comparatively large pion 
mass of M w = 200 MeV significant deviations from the infinite-volume scaling appear only 
for the smallest volume with L = 2 fm. For a physical pion mass of M n = 138 MeV, these 
deviations appear below 4 fm, and for M n = 75 MeV for volume sizes below L = 6 to 8 fm. 

This behavior, however, is very different from the naive expectation for the finite- volume 
behavior: The susceptibility increases in these volumes and is larger than in the infinite- 
volume limit, which cannot be explained in terms of a simple cutoff effect for the long-range 
fluctuations. 

We interpret this behavior as a quark effect in the chiral quark-meson model. This 
behavior is specific to the choice of periodic boundary conditions for the quark fields in the 
spatial directions of the finite volume. As observed in a s yst ematic study of the effects of 
the quark boundary conditions in the quark- meson model 60|], zero-mode effects in a finite 
volume lead to an increase in the chiral quark condensate and a corresponding decrease of 
the pion mass in an intermediate volume size. This affects in turn also the longitudinal 
susceptibility \a and leads to larger values than in the infinite- volume limit. In contrast, 
with a choice of anti-periodic boundary conditions for the quark fields, the lowest momentum 
mode acts as a mass gap, which increases with decreasing volume, and these effects are 
absent (60]. 

It is likely that these quark effects are stronger in the model calculation than in actual 
QCD, where quarks are subject to confinement at low momentum scales. Nonetheless, 
evidence for this effect in QCD comes from lattice QCD simulations and the approach via 
Dyson- S chwinger Equations: A decrease of the pion mass in a finite volume was observed 
in quenched [90J as well as unquenched [9l| lattice simulations, and in a Dyson-Schwinger 
approach to QCD 92|. The authors of the last reference interpret their results as an effect 
of quenching, but cannot exclude that it is also present if all pion effects are taken into 
account. 

Our results demonstrate that a finite-size scaling analysis is feasible for sufficiently small 
volumes. On the other hand, additional quark effects already manifest themselves already 
for larger volume sizes and affect the analysis in terms of infinite- volume scaling analysis, if 
they are not taken into account. 
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Figure 6: Rescaled cr-susceptibility for M n = 138 MeV (left panel) and M w = 200 MeV (right panel) 
as a function of z for various box lengths L. 



V. CONCLUSIONS 



In this paper we have performed a scaling analysis of the chiral order parameter in infinite 
and finite volume with the aid of an effective low-energy model, namely the chiral quark- 
meson model. Using functional RG methods we have computed the scaling functions in 
infinite and finite volume and have analyzed the behavior of the chiral susceptibility in 
detail. In contrast to earlier studies of the scaling behavior of three-dimensional theories, we 
have performed the scaling analysis in d — 4 Euclidean dimensions, where the temperature in 
a field-theoretical sense is determined by the extent of the Euclidean volume in (imaginary) 
time direction. With a such a setup, scaling of the order parameter associated with a three- 
dimensional universal class is only observed when the conditions for dimensional reduction 
are met. We have confirmed that the observation of scaling requires therefore that T C L ^> 1, 
M^L 3> 1 and M n /T c <C 1, where T c is the critical temperature for M w — > and L — > oo. 

In our study of scaling in infinite volume we have only observed scaling behavior of the 
chiral susceptibility and the chiral order parameter for rather small pion masses, M n < 
1 MeV. Corrections to scaling become soon apparent for M n > 1 MeV. Moreover, we have 
compared our result for the scaling function of the chiral order parameter for M n — > with 
the rescaled order parameter for different pion masses with M n > 75 MeV. We have found 
that the results for the rescaled order parameter for these pion masses fall almost on one 
line. This appears to indicate the proper scaling behavior in this pion mass regime but 
only if we disregard the results for M n 1 MeV. This observation might be of interest 
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for scaling studies in lattice QCD simulations, even though non-universal quantities such 
as the normalization constants T and H in our scaling analysis differ from those in full 
QCD. Turning the argument around, this observation and, in consequence, the too-strong 
dependence of the critical temperature on the pion mass might also be considered as a 
warning for studies of the phase-diagram based on (P)NJL/(P)QM-type models. As we 
have argued, this discrepancy can possibly be traced back to the unknown dependence of 
the (UV) parameters on the current quark mass in these models. One way to lift this 
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ambiguity in the parameters is the use of QCD RG flows 

In finite volumes we have also found that the chiral order parameter shows the expected 
scaling behavior. In addition, we have studied the finite-volume effects which appear if a 
conventional infinite- volume scaling analysis is performed in finite volumes with box lengths 
L = 2, . . . , 30 fm and with pion masses M n > 48 MeV. This allows us to depict clearly the 
deviations from infinite-volume scaling behavior which are due to the presence of a finite 
volume. We have found that the chiral susceptibility scales as Xo- ~ L 2 for small volumes, 
M n L < 1. If this condition is not met, then the susceptibility does not scale as Xa ~ L 2 . On 
the contrary, we have shown that the height of the peak of the susceptibility increases in this 
pion mass regime with decreasing volume size. We would like to stress that this is a quark 
effect which might be less pronounced in lattice QCD simulations, but which is presumably 



still present. It has indeed been found in lattice simulations (90|, DSE studies 92( and in an 
earlier RG study 46] that the pion mass exhibits a minimum as a function of the volume size 
depending on the actual value of the dimensionless quantity TL. This effect can be traced 
back to the spatial zero modes of the fermions and it is therefore only present when periodic 
boundary conditions for the fermions are applied 46]. In any case, this volume dependence 
of the pion mass clearly affects the susceptibility. 

Overall, we believe that our results provide further insight into the various aspects of the 
scaling behavior of the chiral order parameter and will provide helpful information for the 
scaling analysis of lattice QCD results. In addition, we think our results will also help to 
further improve the construction of QCD low-energy models. 
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